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ABSTRACT 


This  paper  is  concerned  with  the  detection  of  signal  waveforms  by  a  dis¬ 
tributed  surveillance  network  comprised  of:  a  collection  of  spatially 

separated  sensors,  and  (2)plocal  signal  processors  collocated  with  the  sen¬ 
sors.  The  local  signal  processors  are  assumed  to  implement  likelihood  ratio 
tests  to  detect  the  presence  or  absence  of  the  signals.  Signal  detections  may 
be  used  for  local  decisionmaking  or  passed  upward  to  a  fusion  center  for  fur¬ 
ther  processing.  In  either  case,  the  local  detection  thresholds  cannot  be 
determined  independently,  but  must  be  determined  jointly  to  optimize  overall 
surveillance  system  performance.  Results  are  presented  concerning  the  nature 
of  this  threshold  computation  for  a  number  of  interesting  cases. 
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SECTION  1 
INTRODUCTION 

Classical  detection  theory  (see,  e.g.,  [l]-[4]*)  has  been  motivated  pri¬ 
marily  by  single  sensor  detection  problems.  Although  the  signal  processing 
solutions  of  classical  detection  theory  are  in  principle  equally  applicable  to 
multiple  sensor  detection  problems,  in  practice  these  solutions  may  require 
the  communication  of  raw  received  signals  from  physically  remote  sensors  to  a 
central  processing  location.  In  many  surveillance  systems,  particularly  mili¬ 
tary  systems,  such  communication  capability  is  unavailable  for  reasons  of 
cost,  reliability,  bandwidth,  survivability,  security,  and  similar  factors. 

In  practice,  the  sensors  and  associated  local  signal  processors  in  dis¬ 
tributed  surveillance  networks  tend  to  be  designed  relatively  independently  of 
one  another,  with  communication  restricted  to  higher  level  signal  characteris¬ 
tics  (e.g.,  reports  of  emitter  detections).  While  such  practice  has  the  vir¬ 
tue  of  simplicity,  there  is  a  potential  loss  in  performance  due  to  considering 
each  sensor  individually,  rather  than  as  an  element  in  an  overall  surveillance 
system. 

Motivated  by  such  considerations,  Tenney  and  Sandell  formulated  and 
solved  a  number  of  distributed  detection  problems  (with  random  variable  obser¬ 
vations)  in  [5).  From  a  technical  point-of-view,  these  problems  were  problems 
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in  team  theory  [6]  or  decentralized  control  theory  [7J.  Subsequently,  a  num¬ 
ber  of  authors  have  generalized  the  formulations  in  [5]  in  various  directions 
[ 8 J — 1 12 J ;  however,  all  have  considered  only  the  case  of  random  variable 
observations . 

In  this  paper,  we  generalize  the  distributed  detection  theory  results  of 
[5J  to  the  case  of  waveform  (i.e.,  stochastic  process)  observations.  Thus  we 
will  have  two  hypotheses  to  test, 

1.  H° :  signal  is  absent,  and 

2.  H1:  signal  is  present. 

Corresponding  to  these  two  hypotheses,  we  have  the  potential  decisions  u^-0 
(H°  is  true)  and  uj-1  (H1  is  true)  made  at  the  i-th  local  sensor.  The  cost  of 
decision  errors  is  measured  by  a  cost  function  .Ku^.u^H)  (for  the  case  of  two 
sensors);  it  is  desired  to  choose  local  decisions  u^  on  the  basis  of  locally 
available  observations  yi(t)  to  minimize  the  expected  value  of  J. 

Note  that  in  the  special  case  of  a  surveillance  system  in  which  the  deci¬ 
sions  Ui  are  transmitted  to  a  central  fusion  center  (Fig.  1-1)  where  a  final 
decision  is  made,  the  cost  J  has  the  special  form 

J(Ul,u2,H)  -  J’(f(u1,u2),H)  (1-1) 

where  J'(u,H)  is  the  global  decision  cost  and  f(ulfu2)  is  the  (fixed)  fusion 
law  (e.g.,  voting)  used  to  determine  u  from  u1  and  u2» 

In  the  sequel,  we  will  assume  that  the  local  signal  processing  imple¬ 
mented  at  the  i-th  sensor  consists  of  the  following  two  steps:* 

*As  will  be  noted  subsequently,  only  in  special  cases  is  it  in  fact  optimal  to 
base  the  detection  decision  on  the  local  likelihood  ratio. 
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1.  Determine  the  local  log-likelihood  ratio  tj;  and 

2.  Implement  the  test* 


> 

<  Ti  .  (1-2) 

H0 

We  will  then  determine  the  optimal  thresholds  T^. 

In  Section  2,  we  consider  the  case  in  which  each  y^(t)  is  either  a  known 
signal  plus  noise  (Hl),  or  just  noise  (H°).  In  the  special  cases  in  which: 

(1)  the  noise  is  uncorrelated  between  sensors,  and  (2)  the  noise  is  correlated 
between  sensors  but  is  white  and  the  signals  are  linearly  dependent,  we  are 
able  to  obtain  sufficient  statistics  for  local  decisionmaking.  The  technique 
involves  expanding  the  observation  processes  into  a  Karhuren-Loeve  series. 
Numerical  results  are  provided  for  the  case  of  detection  of  a  sinusoidal  sig¬ 
nal  in  white  Gaussian  noise.  Numerical  results  are  presented  for  both  the 
general  Bayesian  case  and  for  the  case  with  a  centralized  fusion  center. 

In  Section  3  we  consider  the  case  in  which  the  signal  is  a  Guassian  sto¬ 
chastic  process.  Results  are  more  difficult  to  obtain  here,  since  the  log- 
likelihood  ratio  is  a  nonlinear  (quadratic)  function  of  the  observations  even 
in  the  centralized  case.  However,  in  the  important  special  case  of  an  ideal 
band  limited  signal  in  white  noise,  we  are  able  to  obtain  a  closed  form 
(albeit  rather  complex)  expression  for  the  joint  probability  density  function 
of  the  local  log-likelihood  ratios.  This  expression  permits  us  to  solve  for 

*The  notation  means  "choose  H*  if  li>Tj,  and  choose  H®  if  ti<Ti." 
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the  optimal  distributed  detection  thresholds.  Numerical 
for  both  the  general  Bayesian  case  and  for  the  case  with 
center. 

Section  4  contains  a  summary,  some  conclusions,  and 
future  research. 


results  are  presented 
a  central  fusion 

suggestions  for 
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SECTION  2 

KNOWN  SIGNAL  IN  NOISE 

In  this  section  we  consider  a  distributed  detection  problem  in  which  the 
signals  are  known  deterministic  waveforms  and  the  noise  processes  are  colored 
Gaussian  processes.  As  stated  earlier,  we  will  restrict  attention  to  the  case 
of  two  detectors  each  of  which  must  make  a  binary  decision  based  on  local 
observations. 

2.1  PROBLEM  FORMULATION  AND  PRELIMINARY  ANALYSIS 

We  assume  that  there  are  two  sensors,  indexed  by  i-1  ,2  and  that  there  are 
two  hypotheses  to  be  tested  based  on  the  sensor  observations.  The  observa¬ 
tions  under  the  two  hypotheses  are  modeled  by  the  equations 


H1:  yi(0  “  yfii  Si(t)  +  ni(t)  T0<t<Tf 

H°:  yj(t)  -  ni(t)  T0«<Tf  . 


(2-1) 


We  assume  that  the  S£(t)  are  known  signals  with  unit  energy  and  are  zero  out¬ 
side  the  interval  [0,T]  where  To<0<T<Tf.  The  ni(t)  are  assumed  to  be  zero- 
mean  Gaussian  processes  where 


Etn^On-jCt)}  -  Kjj(t,t)  ,  i,j«l,2  (2-2) 

and  we  assume  (to  avoid  the  possibility  of  singular  detection) 
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Kii(t.T)  -  iq(t,T)  +  Nifi(t-t)  (2-3) 

with  Nj^O.  Note  that  we  do  not  assume  that  the  noise  processes  are  uncorre- 

n 

lated  between  sensors,  i.e.,  we  allow  K^j(t,T^  to  be  nonzero  for  i*j. 

The  approach  we  will  use  to  determine  the  optimal  distributed  detection 
laws  consists  of  expanding  the  received  waveforms  in  a  Karhunen-Loeve  (K-L) 
expansion  [1]  and  considering  the  problem  formed  by  truncating  the  infinite 
series  of  K-L  coefficients  to  the  first  K  terms.  The  truncated  problem  can  be 
approached  via  the  technique  of  [ 5 J  and  results  for  the  waveform  problem 
obtained  by  taking  the  limit  as  K-mo. 

Thus  we  expand  yf(t)  via* 


K  k  k 

yi(t)  -  lim  E  yi$i(t)  .  T0<t<Tf 
K->«  k-1 


k 

where  the  <J>i(t)  satisfy 


k  k  "^f  c  k 

Xi<t>i(t)  =  /  Ki(t  ,u)if>i(u)du  ,  T0<t<Tf  . 

To 


Under  H*  we  have 


(2-4) 


(2-5) 


k  k  k 

yi  -  si  +  n^  (2-6) 


where 


*x*lim  Xfc  is  defined  to  mean  lim  lx-xyll2“0. 
k>«°  k+« 
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k  Tf  k 

si  ^  I  Si(t)^(t)dt 


(2-7) 


k  Tf 

ni  A  J  nl(t>4,i(t)dt  , 

Tn 


(2-8) 


while  under  H°  we  have  that 


It  k 
yt  -  ni 


(2-9) 


It  is  straightforward  to  show  that 


E{n^}  *  0 


(2-10) 


k  2 

E{(n1)2)  -  Ai  +  % 


(2-1 L ) 


and  that  the  K-L  coefficients  corresponding  to  each  sensor  are  uncorrelated: 


Z  k 

E{nini)  *  0  £*k 


(2-12) 


However,  the  K-L  coefficients  for  the  two  sensors  are  correlated,  that  is, 


Tf  Tf 

E{(yJ-E{y^|Hj})(y|-E{y^|Hj}|HJ}  *  /  /  4>^(t)K^2(t ,u) ^(u)dtdu  £  c*k  , 


To  T0 


j-0,1  ,  i*k 


(2-13) 


which  is  not  zero  in  general. 
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Using  the  results  of  [5j  we  can  determine  (implicitly)  the  optimal  dis¬ 
tributed  decision  law  for  the  problem  based  on  the  first  K  coefficients  of  the 
K-L  expansion.  Define 


K  12  K 

2i  *  lyi.yi.*** .y i 1 


(2-14) 


K  K 

Ai(yi)  =  n  exp 
k=l 


k  5  k  k 
(Si)  -2yiSi 

2(Xk+Nj) 


(2-15) 


The  optimal  decision  law  is  then  given  by  the  solution  of  the  following  two 


equations : 


V4}  < 


y^ju  =0 
2  2 


<  aj3+aj4  /.../  p(y^|y^H°)dZK 


AJ5+AJ2  /.../  pCz^lzl.Hild^ 

yK|u  *0 

0  zi'  i 

A  (yK)  2  •  — •  -  -  -  —  *  *  1 — " — — — 

2  2  <  AJg+AJl4  /.../  P<^|y|,H°)d^ 


(2-lba) 


(2-16b) 


where  the  AJn  are  given  by  the  following  equations: 
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AJ1  -  JCO.l.H1)  -  JU.l.H1)  (2-17a) 

AJ2  =  JCO.O.H1)  +  JCl.l.H1)  -  Jd.O.H1)  -  J(0,l,Hl)  (2-17b) 

AJ3  *=  J(  1 ,1  ,H°)  -  J(0,1  ,H°)  (2-17c) 

AJ^  >=  J(1,0,H°)  +  J(0 ,1  ,H°)  -  J(O,O,H0)  -  J(1,1»H°)  (2-17d) 

AJg  -  Jd.O.H1)  -  Jd.l.H1)  (2-17e) 


and 

AJg  =  J(1,1,H°)  -  Jd.O.H0)  .  ( 2— 1 7 f ) 

Note  that  the  decision  law  for  each  sensor  is  required  to  determine  the 
region  of  integration  for  the  right-hand  side  of  Eq.  2-16.  Thus,  the  deter¬ 
mination  of  the  optimal  distributed  decision  law  requires  the  solution  of 
coupled  nonlinear  functional  equations.  No  general  analytic  solution  to  these 
equations  has  been  determined  and  numerical  techniques  do  not  seem  to  be  com¬ 
putationally  feasible.  In  general,  it  is  necessary  to  assume  some  special 
form  for  the  local  signal  processing  (e.g.,  a  likelihood  ratio  computation) 
and  optimize  the  thresholds.  However,  in  certain  special  cases  more  can  be 
said. 
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2.2  SPECIAL  CASES 
2.2.1  Uncorrelated  Noise 

The  easiest  case  to  consider  is  that  in  which  the  noise  is  uncorrelated 
between  sensors,  i.e.,  Kj2(t,r)=0.  Under  this  assumption  we  have  that 


K  K  K 

P(lil2j*H*)  =  P<Xi|H*)  i  ,j«l  ,2,i*j  ,1=0,1 
and  thus  Eq.  2-16  becomes 


VZi> 


AJj+AJ^  /.../  p(z£)Hl)d2£ 


U  2  2 

>  -  A  TK 

<  A J  ^+A J  ^  J...J  p(x||H°)d^K  1 

^2  I u2“° 


(2-18) 


(2-19) 


AJ5+AJ2  J...J  P(ljlHl)dij 


A2(yK)  > 


£K|u  =0 
I  1 


<  AJ3+AJv  /.../  p(yK|H°)dZK 

Z^iui“0 


A 

—  2 


(2-20) 


Note  that  the  right-hand  sides  of  the  above  equations  do  not  depend  upon 

K 

the  local  observations,  but  are  constants,  which  we  define  as  Tj.  Thus,  the 
optimal  distributed  decision  law  is  given  by  a  pair  of  local  likelihood  ratio 
tests  where  the  thresholds  are  implicitly  defined  by  Eqs.  2— 19  and  2-20. 

If  we  take  the  limit  as  then  one  can  show  (see  (lj)  that  the  optimal 
decision  laws  can  be  written  as 
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I  yi(t)gi(t)dt  >  -  /  Si(t)gi(t)dt  -  Ti  (2-21) 

<  2  *  T 


where  the  left-hand  side  of  Eq.  2-21  is  determined  by 


gi(t)  -  /  Qi(t ,u)si(u)du 


(2-22) 


and  Qi(t,u)  is  defined  by 


/  Kii(t ,u)Qi(u,v)du  -  6( t-v ) 


(2-23) 


The  right-hand  side  of  Eq.  2-21  is  a  constant  (observation-independent) 
threshold  and  the  Ti  satisfy  the  following  nonlinear  coupled  algebraic 
equations : 


t1  *  in 


(  AJ  +AJ  [  1-erf  ( (  t  +m  )/^2ra  ) )  ) 

)  1  2 _ 22  2  ( 

|  AJ3+AJ4ll-erf((T2-m2)/^2m2)J  j  p° 


(2-24) 


t2  *  In 


(  AJ  +AJ  [1-erf ((t  +m  )/\/2m  )J  ) 
15  2  1  1  w  1  |  pi 

|  AJ6+AJH[l-erf((t1-m1)/y2m1))  j  p° 


(2-23) 


where 


Vet  tf 

mi  A  -  /  Si(t)gi(t)dt  . 


2  Tn 


(2-26) 
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Note  that  in  general,  Eqs.  2-24  and  2-25  do  not  have  a  unique  sol,  tion 
and  thus  the  optimal  decision  law  must  be  determined  by  computing  all  solu¬ 
tions  and  then  comparing  their  performance.  The  decision  law  of  Eq.  2-21  ex¬ 
tends  the  results  of  [5]  to  the  case  of  waveform  observations  and  this  law  can 
be  interpreted  as  a  pair  of  local  likelihood  ratio  tests  with  thresholds  which 
are  jointly  optimized  according  to  a  system-wide  measure  of  performance. 


2.2.2  Linearly  Dependent  Signals  and  Correlated  White  Noise 

In  this  special  case  we  assume  that  s ,“S-As  so  that  the  observations  are 
given  by 


H1:  yi(t)  -  s(t)  +  nj (t)  0<t<T 

H° :  y±(t)  ■=  nt(t)  0<t<T  . 


(2-27) 


Furthermore  we  assume  that  the  n^Ct)  are  zero-mean  unit  spectral  height,  white 
Gaussian  noise  processes  where 


E{n1(t)n2(r)}  -  p6(t-t) 


(2-28) 


As  usual  we  expand  the  received  waveform  in  a  K-L  expansion  where  now  we 
k 

choose  the  4>i(t)  to  be  any  complete  orthonormal  set  such  that 


^  m  4>2(t)  “  *k(t)  0«<T 


(2-29) 


and 


♦  l(t)  -  s ( t )  0<t<T  . 

It  is  easy  to  verify  that  if 


(2-30) 
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K  k  k 

yi(t)  -  lim  l  y i 4>  (t) 
K+«  k«l 


(2-31) 


then 


1  (  +  wi  :  h1 

^  "  |  Wi 


(2-32) 


k  |  :  H» 

yi  *  ;  k*2,...  (2-33) 

I  Vi  :  H° 


where 


T 

wi  A  /  <t>i(t)ni(t)dt 

”  0 


and 


E{wkw|}  »  p6(k-i)  .  (2-34) 

k 

Note  that  the  coefficients  yi  for  k>2  have  the  same  statistics  under  both 
hypotheses  and  are  independent  of  the  y|.  Thus  any  optimal  decision  law  will 
use  only  the  y1  to  decide  between  the  hypotheses  H®  and  H1 .  The  y^  are  thus 
scalar  sufficient  statistics  for  optimal  decisionmaking  in  the  team  decision 
problem  with  information  pattern  defined  by  Eqs.  2-27  and  2-28. 

Any  optimal  decision  law  can  therefore  be  specified  for  this  problem  by 
defining  the  regions  of  the  real  line  in  which  y*  must  lie  for  a  decision  of 
H1  to  be  made.  These  regions  can  be  specified  by  their  endpoints,  and  thus 
the  decision  law  can  be  characterized  by  a  set  of  endpoints,  or  thresholds, 
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for  each  detector.  The  optimal  distributed  decision  law  is  determined  by 

24-1  24  Li 

choosing  collections  of  thresholds  {[Ti  ,Ti  i=l,2  as  follows: 

21  ,T  24+1 

(  1  if  Tj  <  /  Si(t)yi(t)dt  <  Ti  ,  4=1 . L 

° 

Ui  =  \  .  (2-35) 

\  0  otherwise 

Necessary  conditions  for  the  optimality  of  the  thresholds  are  readily 
developed  if  Li  is  specified.  Since  the  number  Li  of  such  thresholds  is  arbi¬ 
trary  one  must  compare  the  performance  of  laws  with  different  Li's  to  deter¬ 
mine  the  optimal  distributed  detection  law.  We  have  calculated  optimal  laws 
based  on  the  assumption  that  Lj*L2=L=1  or  3.*  In  all  cases  the  laws  based  on 
the  assumption  that  L«1  proved  superior.  While  we  have  been  unable  to  prove 
that  L=1  is  optimal,  we  will  only  consider  that  case  in  the  sequel. t 
The  optimal  detection  law  for  the  L*1  case  is  given  by 

T  1 

/  s(t)yi(t)dt  >  Ti  (2-36) 

0  < 

0 


wherre  Tj  and  T2  satisfy 


♦Clearly  ui»0  as  y|-»— «°  and  ui=l  as  y|-*+-»  so  that  an  even  number  of  intervals 
is  required  and  thus  L  must  be  odd. 

tNote  that  the  necessary  conditions  for  optimality  consist  of  2L  nonlinear 
coupled  equations;  the  assumption  of  L=1  leads  to  considerable  reduction  in 
the  computational  burden. 
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/eT  (  aj 

VJ  ±_  !  )  1 

2  "j 


+AJ  erf 
2 


UY>(Yp(Y>/V)/''^  J 

I 


(2-37) 


AJ3+AJ4  erf[(T2-pT1  )/\l-p2l 


T 


2 


AJg+AJ2  erf  [(Ti-<^F'-p(T2-^r’))/Vl-P2  j 

2  j 

AJg+AJ4  erf  [(T1-pT2)/Vl-P2)  j 

For  the  special  cases  considered  above  the  "local"  likelihood  ratio  is  a 
sufficient  statistic  for  optimal  distributed  decisionmaking.  More  generally, 
coupled  functional  equations  must  be  solved  to  determine  the  optimal  distrib¬ 
uted  detection  law.  However,  given  the  difficulty  of  solving  such  equations 
and  the  above  results,  a  reasonable  choice  for  a  distributed  detection  law  is 
to  use  local  likelihood  ratio  tests  with  globally  optimized  thresholds.  Note 
that  this  choice  implies  that  the  same  signal  processing  can  be  used  whether  a 
sensor  is  used  alone  or  as  part  of  a  surveillance  network  -  a  strong  practical 
requirement!  In  Section  3  we  will  take  this  approach  when  considering  the 
detection  of  an  unknown  signal  in  noise. 

2.3  NUMERICAL  EXAMPLES 

Here  we  consider  an  example  involving  a  scalar  observation  with  corre¬ 
lated  noise.*  First  we  introduce  the  observation  model  and  rewrite  the  neces¬ 
sary  conditions,  then  in  subsections  2.3.1  and  2.3.2  we  consider  two  different 
global  cost  functions. 


*The  special  case  of  uncorrelated  noise  reduces  to  the  scalar  observation 
problem  studied  in  [5]  and  thus  we  do  not  consider  it  here. 
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Assume  that  the  observations  are  given  by 

H1 :  Yi(t)  ■  ^2E^  sin(2»t)  +  n^(t)  0<t<l 

H°:  yi(t)  *  nj(t)  0<t<l 

where  n^(t)  is  zero-mean  unit  variance  white  Gaussian  noise  with 

E{n1(t)n2(r)}  -  pfi(t-T)  . 


(2-39) 


(2-40) 


We  assume  p0*pl=l/2  and  note  that  the  detection  law  can  thus  be  written  as 

i  -  T  1 

y  A  2  /  yi(t)sin(2irt)dt  >  (2-41) 

1  0  < 

0 

where  the  Tj  are  solutions  to  Eqs.  2-37  and  2-38. 

2.3.1  Bayesian  Formulation 

We  assume  that  the  thresholds  are  to  be  selected  so  as  to  minimize  the 
expected  Bayes  cost  where 


/  0  if  u  *u2«H 

J  -  !  1  if  u^u2  .  (2-42) 

\  k  if  u^u-^H 
1  2 


This  type  of  cost  criterion  arises  in  situations  where  it  is  not  precisely 
twice  as  costly  to  make  two  errors  as  it  is  to  make  one.  For  example,  if 
weapons  are  collocated  with  sensors  and  targets  are  automatically  attacked 
when  detected,  then  having  two  missed  detections  may  be  much  more  serious  than 
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having  one  missed  detection.  This  situation  would  be  modeled  by  choosing  k>2, 
so  that  a  double  error  is  more  than  twice  as  costly  as  two  single  errors. 

For  J  defined  by  Eq.  2-42,  Eqs.  2-37  and  2-38  become 


T 


1 


I  l+(k-2)erf 


in< 


t -JT) 
2  V  2  1  V  l 

^  1-p2 


(k-1 )-(k-2)erf 


T  -pT 
2  1 


(2-43) 


(2-44) 


and  we  note  that  the  thresholds  for  the  locally  optimal  test*  (T VV2' 
satisfy  Eqs.  2-43  and  2-44.  This  can  be  seen  by  noting  that 


*That  is  the  test  with  thresholds  that  would  be  optimal  for  each  sensor  con¬ 
sidered  in  isolation  when  a  minimum  probability  of  error  criterion  is  used. 


18 


ALPHATECH,  INC. 


-  ln(l)  -  0  .  (2-45) 


A  a 

Equations  2-43  and  2-44  are  symmetric;  if  CT^.T^)  is  a  solution,  then 

A  A 

(^E^-Tj,  is  als0  3  solution’  This  is  seen  by  notin8  that  if 


then 
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Vv1 


, 


v5: 


in 


(k-l)-(k-2)erf 

A  A 

(t  -pT  ')/\ll- p2 

2  1  * 

|n-(k-2)erf 

_<vV= 

.2-p(T1-^))/^l-p2 

^  , 


Vs 


in< 


( k— 1 )-(k-2) 

1-erf 

l+(k-2) 

1-erf 

<  Vvt2)-p(Vvti  y 

vs 


in , 


^l+(k-2)erf 

A  A 

<>/yv>/y»(Vvy-Vv 

^^2 

)  <k-l)-(k-2)erf 

(  Je~-T  )-p(  /r-T  ) 

\  2  2  V  1  1 

(2-47) 


The  costs  associated  with  these  two  symmetric  solutions  are  identical  and 
thus,  if  a  solution  with  Ti*^Ei/2  exists,  only  one  of  this  pair  need  be 
evaluated . 

Graphical  solution  of  Eqs.  2-37  and  2-38  shows  that  there  are  at  most 
three  solutions  to  these  equations,  and  for  certain  values  of  k  and  p  there  is 
only  one  solution  (the  "locally  optimal  test"  (LOT)).  In  Figs.  2-1  and  2-2  we 
plot  the  Bayes  cost  associated  with  the  decentralized  likelihood  ratio  test 
(DLRT)  and  with  the  LOT  solutions  for  ^E^“l ,  ^Ej«2  and  various  p  and  k. 
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Figure  2-1  is  a  plot  of  the  Bayes  cost  as  a  function  of  k  for  p*>0  and 
p«0.5.  Note  that,  for  larger  p,  the  optimal  DLRT  solution  can  be  much  better 
than  the  LOT.  This  is  because  the  optimal  solution  skews  the  thresholds  to 
avoid  the  double  errors  which  are  more  common  for  larger  p  and  more  costly  for 
larger  k. 

Note  that  for  Che  LOT  the  cost  increases  as  an  affine  function  of  k  since 
the  probability  of  a  double  error  is  independent  of  k.  The  expected  cost  for 
the  DLRT  never  exceeds  1  since  that  value  can  be  obtained  by  a  suboptimal  law 
with  thresholds  set  so  that  one  detector  always  decides  0  while  the  other 
always  decides  1. 

Figure  2-2  is  a  plot  of  the  Bayes  cost  as  a  function  of  p  for  k*5.  Again 
we  see  that  as  p  increases  the  DLRT  becomes  much  better  than  the  LOT.  This 
occurs  because,  as  p-*T ,  the  probability  of  a  double  error  (for  fixed  thresh¬ 
olds)  increases.  The  DLRT  solution  skews  the  thresholds  to  decrease  the  prob¬ 
ability  of  a  double  error.  This  yields  a  27  percent  decrease  in  cost  over  the 
LOT. 

These  results  indicate  that  in  some  cases  there  is  a  significant  gain  to 
be  had  by  using  the  optimal  decentralized  likelihood  ratio  test  rather  than  a 
naive  approach  which  ignores  the  correlation  between  sensors. 

2,3.2  Surveillance  System  Design 

In  this  subsection  we  assume  that  a  local  decision  u^  based  only  on  the 
information  provided  by  sensor  i  is  sent  to  a  fusion  center  where  a  global 
decision  u  is  made.  We  assume  that  the  fusion  center  decision  rule  is  defined 
as  u-1  if  and  only  if  u^u^-l.  We  display  the  overall  performance  of  such  a 
surveillance  system  via  a  generalized  network  receiver  operating 
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characteristic  (ROC)  curve,  which  is  a  plot  of  the  surveillance  network's 
probability  of  detection  versus  probability  of  false  alarm.  In  addition  to 
plotting  the  performance  of  the  DLRT  we  also  plot  the  performance  of  the  opti¬ 
mal  centralized  detection  law.  This  allows  us  to  determine  how  the  perfor¬ 
mance  of  the  surveillance  system  is  degraded  by  requiring  that  only  local 
decisions  (processed  information)  rather  than  raw  sensor  data  be  transmitted 
to  the  fusion  center. 

The  ROC  can  be  obtained  by  varying  the  ratio  of  the  cost  of  a  false  alarm 
to  the  cost  of  a  missed  detection.  If  we  let  the  false  alarm  cost  be  unity 
and  the  missed  detection  cost  be  a  then  the  necessary  conditions  for  optimal¬ 
ity  become 


.£.j_ 

T 

j  1-erf  ^(Ti-^r-p(T2->|r))/VT:V 

2  ft 

|  1-erf  [jTj-pT^/V^P1^ 

■j*2  1  ( 

T  n  ■  _  ____  On  ' 

1  n  -  -  -  xn  • 

2  ^  { 

1-erf  Q^-pTjJ/VlVj 

Figure  2-3  depicts  the  performance  of  the  optimal  centralized  test  and 
that  of  the  DLRT  for  the  case  where  the  sensors  are  identical, 

We  see  that  as  p  increases  the  centralized  and  DLRT  results  become  more  and 
more  similar.  Figure  2-3  also  shows  that  as  p  increases  the  performance 
degrades. 
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Note,  however,  that  the  performance  of  the  DLRT,  while  more  closely 
approaching  that  of  the  centralized  test,  degrades  as  p+1 .  This  can  be  under¬ 
stood  by  noting  that  if  p  is  the  probability  of  a  local  decision  being  wrong 
then  the  probability  of  a  double  error  is  approximately  p2  when  p«0  but  is  p 
when  p*l.  Since  the  fusion  center  always  makes  an  incorrect  decision  when 
both  local  decisions  are  wrong,  the  performance  degrades  as  p-*-l. 

The  phenomenon  of  the  DLRT  and  centralized  results  growing  closer  togeth¬ 
er  as  p  increases  is  not  universal.  Figure  2-4  illustrates  the  behavior  of 
these  two  decision  laws  for  the  case  of  asymmetric  sensors,  and  Vv2- 

We  see  that  for  the  DLRT,  the  performance  degrades  as  p+1 .  The  reason  is 
exactly  the  same  as  in  the  case  of  VvVv  1.  For  the  centralized  case  how¬ 
ever,  the  performance  becomes  perfect  as  p+1 .  This  occurs  because  by  differ¬ 
encing  the  sensor  observations  one  has 

H1:  Ay(t)  -  yx(t)  -  y2(t)  «  (  y/2-1  )sin(2irt)  +  n^t)  -  n2(t)  (2-50) 

H° :  Ay(t)  «=  y^t)  -  y2(t)  «  n^t)  -  n2(t)  .  (2-51) 

As  p+1 ,  Ay+(Y/2-l)sin(2TTt)  if  H1  is  true  and  Ay+0  if  H°  is  true,  thus 
perfect  detection  is  possible. 

These  graphs  illustrate  that  the  performance  difference  between  a  DLRT 
and  a  centralized  test  is  strongly  dependent  on  the  problem  being  considered. 
Similarly,  the  benefit  in  reduction  of  communications  requirements  achieved  by 
using  a  DLRT  rather  than  a  centralized  test  is  highly  problem  dependent.  It 
is  thus  the  case  that  detailed  analysis  is  required  to  determine  whether  a 
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SECTION  3 

UNKNOWN  SIGNAL  IN  NOISE 


In  this  section  we  consider  a  distributed  detection  problem  in  which  the 
signal  is  a  random  Gaussian  process  and  the  noises  are  white  Gaussian  proces¬ 
ses.  Again  we  consider  only  the  case  of  two  detectors  and  binary  hypotheses. 

We  assume  that  the  i-th  sensor's  observation  under  the  two  hypotheses  is 
modeled  by: 


H1:  yi(t)  -  CiSjXt)  +  n^(t)  0<t<T 

H°:  yi(t)  =  nj  ( t )  0<t<T 


(3-1) 


where  the  n^(t)  are  independent  zero-mean  unit  spectral  height,  white  Gaussian 
noise  processes  and  the  signals  s*(t)  are  zero-mean  unit  power  Gaussian  pro¬ 
cesses  with  known  covariances: 


E{si(t)sj(t)}  &  Kij(t,x)  CKt,t<T  .  (3-2) 

s 

If  we  expand  the  yi(t)  in  K-L  series  (generated  by  Kii(t,t))  then  the  op¬ 
timal  distributed  detection  test  is  given  by  Eqs.  2-16  and  2-17,  where  the 
A±(y±)  are  now  quadratic  functions  of  the  observations  yj(t).  Unlike  the 
known  signal  in  noise  case,  all  the  K-L  coefficients  have  statistics  which 


28 


ALPHATECH,  INC 


depend  upon  which  hypothesis  is  true.  Thus  none  of  the  simplifications  possi¬ 
ble  in  subsection  2.2  can  be  applied  -  to  determine  the  optimal  detection  law, 
coupled  nonlinear  functional  equations  must  be  solved. 

Since  we  cannot  solve  the  equations  defining  the  optimal  distributed 
detection  test  we  look,  instead  for  a  good  alternative.  The  test  we  consider 
is  that  motivated  by  the  previous  section:  the  distributed  likelihood  ratio 
test  (DLRT).  These  tests  are  defined  by  the  equation 

1 

*i  >  Ta  ,  i-1,2  (3-3) 

0 


where  is  the  (local)  log-likelihood  ratio  and  the  Tj  are  optimized  for  the 
best  global  surveillance  system  performance.  This  class  of  tests  is  not  easy 
to  analyze  for  the  problem  of  detecting  an  unknown  signal  in  noise.  The  dif¬ 
ficulty  arises  because  the  local  log-likelihood  ratios  are  not  Gaussian.  In 
the  discussion  below  we  consider  a  problem  for  which  we  are  able  to  obtain 
results . 


3.1  IDEAL  BANDLIMITED  SIGNAL 

Consider  a  problem  in  which  the  observations  are  modeled  by* 
H1 :  yi(t)  -  Cis(t)  +  nj:(t)  0<t<T 
H° :  yi(t)  -  ni(t)  0<t<T 


(3-4) 


*The  bandpass  version  of  this  problem  is  readily  treated  in  an  entirely  analo¬ 
gous  fashion;  such  a  formulation  is  a  reasonable  model  for  passive  detection 
of  a  radio  transmitter  operating  at  a  known  frequency  and  with  a  known 
bandwidth. 
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where  s(t)  is  a  zero-mean  Gaussian  stochastic  process  of  unit  power  with  an 
ideal  bandlimited  spectrum  (Fig.  3-1)  and  the  n^(t)  are  zero-mean  unit  spec¬ 
tral  height,  Gaussian  white  noise  processes.  We  denote  the  covariance  of  the 
signal  by  Ks(t,t). 


S(f) 


R-094C. 


Figure  3-1.  Signal  Spectrum. 

To  determine  the  performance  of  the  DLRT  we  first  determine  the  form  of 
the  local  likelihood  ratios  and  then  calculate  their  joint  probability  distri¬ 
bution  function  (conditioned  on  each  hypothesis). 

3.1.1  Local  Log-Likelihood  Ratio 

If  we  expand  the  received  signal  via  a  K-L  series, 

K 

yi(t)  *  lim  l  y^^Ct)  0<t<T  (3-5) 

K+®  k-1  i 


where  the  $k(t)  satisfy 
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T 

Xk$k(t)  -  /  Ks(t ,u)$k(u)du  0<t<T  (3-6) 

0 


it  is  easy  to  show  [2]  that  the  log-likelihood  ratio  based  on  the  first  K 
coefficients  at  the  i-th  sensor  is  given  by 


i 


1 

2 


K 

l 

k-1 


l+c2Xk 
i 


(yk)2 

i 


1  K 

Z  in(l+c2Xk) 
2  k=l  i 


(3-7) 


Since  the  thresholds  are  constants  in  Eq.  3-3  we  shall  henceforth  work  with 

only  the  data  dependent  portion  of  iK,  which  we  define  as 

1 


A* 

i 


K 

Z 

k-1 


2,k 
c  X 


l+c2Xk  J 
i 


(yk)2 


(3-8) 


In  the  limit  as  K+®  we  obtain  (21 


1  T  T 

*i  *  "  /  /  yj(t)hi(t ,u)yi(u)dudt  (3-9) 

2  0  0 


where  h^(t,u)  satisfies 


r  2  7 

t»i(t,u)  +  /  c  hi(t,z)Ks(z,u)dz  -  c  Ks(t,u)  ,  0<t ,u<T  .  (3-10) 

0  1  i 
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3.1.2  Joint  Probability  Distribution  Function 

The  derivation  of  expressions  for  the  joint  probability  distribution 
functions  of  ij  and  i conditioned  on  H°  and  H1  is  complex  and  thus  the 
details  are  relegated  to  Appendix  A.  Here  we  briefly  sketch  the  derivation. 

We  first  compute  the  joint  moment  generating  function  of  the  iK.  This  is 
relatively  straightforward  as  the  y^  are  (for  each  i)  conditionally  indepen¬ 
dent.  We  then  take  the  limit  as  K-m».  If  we  define 


-rl  -si 

Pi(r,s)  A  E {e  | H^- } 


(3-11) 


then  we  have 


where 


in  u0(r,s) 


-  E  in 
2  k=l 


,  ,  2.k 
1+c  X 


l+(l+r)c2xkj 


2  k  \ 
1+c  XK 


l+(l+s)C2x^ 


(3-12) 


in  Uj (r  ,s) 


-  E  in 
2  k-1 


(1-P2) 


(  l+r(  l-p2)c2x'c)(  l+s(  l-p2)c2X*t)-p2 
k  1  k  2  V 


(3-13) 


c  c  X 
1  2 


Pk  A 


V<l+cfXk)(l+c|Xk) 


(3-14) 
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If  the  observation  time  is  long  in  comparison  to  the  6ignal  time  con¬ 
stants  (i.e.,  if  WT  is  large  enough)  we  have  [1J  that  the  infinite  sums  in 
Eqs.  3-12  and  3-13  can  be  replaced  with  integrals  in  which  the  eigenvalue  mag¬ 
nitude  appearing  in  Eqs.  3-12  and  3-13  is  replaced  by  the  signal's  spectral 
height.  This  yields  (for  the  ideal  bandlimited  signal)  a  trivial  integral  and 
we  have 


l+c2/2W 

WT 

l+c2/2W 

1 

2 

l+( l+r)c2/2W 

l+(l+s)c|/2W 

(3-15) 


Pj(r,s) 


(l+r(l-p2)c2/2W)(l+s(l-p2)c2/2W)-p2 


(3-16) 


where 


c  c  /2W 
1  2 

P  *  —  ~  . .  .  .  (3-17) 

V(1+c2/2W)(1+c2/2W) 


We  can  then  invert  each  moment  generating  function  to  obtain 


p(i1,»2|H°) 


(l  l  ) 

1  2 ' 


Lt 2-1  _ 


(ru/2))2 
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2 


(3-18) 
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[d-P2>Va] 


-A/  2  - 


l  l 

1  2 


(i-p2)y1  U-p2)y2 
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fi  l  (1-p2)y  y 
\  1  2  12 
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\ 
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l 

(l-p2)Y1Y2 

\  1  2  1 

(3-19) 


where  we  have  defined  A=2WT  as  the  observation  time-signal  bandwidth  product, 
Ei«c2T  as  the  signal  energy  received  by  the  i-th  sensor  and  Yi“El/A  as  the 
signal-to-noise  ratio  in  the  signal  bandwidth  (recall  the  noise  had  unit  spec¬ 
tral  height)  and  where  Iv(*)  is  the  modified  Bessel  function  of  order  v. 

Surprisingly,  both  Eqs.  3-18  and  3-19  can  be  integrated  analytically  when 
A/2-1  is  a  nonnegative  integer.  We  thus  obtain  under  this  assumption  the 
joint  distribution  functions: 
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where 


Ti 

Bi  ,  -  . 

( l-p2)c2 


(3-23) 


Via  Eqs.  3-20  and  3-22  we  can  evaluate  the  performance  of  the  DLRT  for 
any  pair  of  thresholds.  It  is  thus  straightforward  to  determine  via  numerical 
techniques  the  optimal  thresholds  for  any  given  problem.  Results  are  given  in 
the  next  section. 

3.2  NU11ERICAL  EXAJ1PLES 

Here  we  will  present  numerical  results  for  the  distributed  detection  of 
an  ideal  bandlimited  signal  in  white  Gaussian  noise.  Both  the  general 
Bayesian  case  and  the  case  of  surveillance  system  design  will  be  considered. 

3.2.1  Bayesian  Formulation 

Using  the  cost  function  defined  in  subsection  2.3.1  we  obtain  the  results 
depicted  in  Figs.  3-2  and  3-3.  In  Fig.  3-2  we  have  plotted  the  expected  cost 
of  the  optimal  DLRT  and  that  of  the  locally  optimal  test*  (LOT)  as  functions 


*For  the  cases  considered  here  the  sensors  are  identical  and  we  define  the  LOT 
as  the  DLRT  satisfying  the  additional  constraint  that  Tj-T^ 
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of  the  double  error  cost  k.  The  different  curves  correspond  to  different 
values  of  signal  energy  (Ei*l,  10,  20,  40,  100),  where  the  time-bandwidth 
product  has  been  held  constant  (A* 10). 

We  note  that  identical  results  are  obtained  for  k<2 ,  but  that  the  costs 
associated  with  the  LOT  become  affine  for  large  enough  k.  This  occurs  because 
when  T1«T2  the  probability  of  a  double  error  cannot  be  reduced  to  zero.  For 
large  enough  k,  the  cost  of  double  errors  dominates  the  growth  of  the  expected 
value  of  the  cost  EJ  and  the  affine  curves  result.  The  cost  for  the  optimal 
DLRT  of  course  never  takes  on  a  value  greater  than  unity  (since  taking  Uj  =  l 
and  u2=0  gives  J<1).  As  the  signal  energy  is  increased  the  performance  of  the 
DLRT  improves  and  the  point  at  which  the  LOT  is  not  longer  globally  optimal  is 
also  increased.  This  latter  effect  occurs  because,  as  the  signal  energy  in¬ 
creases,  the  probability  of  a  double  error  decreases  rapidly.  This  means  that 
the  cost  is  dominated  by  the  cost  of  single  errors  until  k  becomes  very  large. 

In  Fig.  3-3  we  again  plot  the  performance  of  the  LOT  and  the  DLRT  versus 
k,  but  now  signal  energy  is  held  constant,  Ei«20,  and  the  time-bandwidth  prod¬ 
uct  is  allowed  to  vary  (A=2,  10,  100).  Here  we  see  the  effect  of  frequency 
diversity  -  for  small  k  the  LOT  is  optimal  and  the  DLRT  performance  for  A=10 
is  better  than  either  A=2  or  A*100  (as  also  in  Figs.  3-4  and  3-5).  As  k  in¬ 
creases,  the  effect  of  double  errors  increases  and  DLRT  performance  is  domi¬ 
nated  by  the  need  to  minimize  double  errors.  For  large  k  the  DLRT  skews  the 
thresholds  sc  that  one  detector  is  likely  to  have  u**0  and  the  other  have  u*l. 
This  strategy  is  so  effective  that,  for  the  cases  of  Fig.  3-3  with  k>20,  the 
cost  of  double  errors  is  less  than  5  percent  of  the  total  cost.  The  structure 
of  the  Bayes  cost  EJ  for  large  k  is  thus  determined  by  the  probability  of 
single  errors  when  the  thresholds  are  skewed  to  avoid  double  errors. 
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3.2.2  Surveillance  System  Design 

Now  we  present  results  for  the  case  described  in  subsection  2.3.2,  i.e., 
the  case  in  which  a  global  decision  is  made  at  a  fusion  center  by  "adding" 
together  the  local  decisions  made  from  the  observations  of  the  individual 
sensors. 

In  Figs.  3-4  and  3-5  we  plot  the  probability  of  a  miss  (i.e.,  that  H°  is 
decided  when  H*  is  true)  as  a  function  of  the  time-bandwidth  product  A*2WT, 
for  fixed  probability  of  false  alarm  (i.e.,  that  H1  is  decided  when  H°  is 
true).  In  Fig.  3-4  we  have  set  Pp=0.1  and  in  Fig.  3-5  we  have  set  Pp=0.001. 
The  various  curves  represent  the  performance  of  the  DLRTs  and  centralized  LRTs 
for  different  signal  energies. 

The  curves  labeled  Ej  plot  the  DLRT  performance  for  two  sensors  each  of 
which,  under  H1 ,  observes  a  signal  of  energy  E<j.  The  curves  labeled  Ec  plot 
the  performance  of  the  optimal  centralized  likelihood  ratio  test  where,  under 
H1,  the  total  signal  energy  received  is  Ec.*  In  both  figures  we  see  the  an¬ 
ticipated  ranking  of  performance  curves:  Ec<E,j<2Ec.  This  occurs  since  the 
DLRT  for  the  two  sensor  case  must  perform  at  least  as  well  as  the  one  sensor 
centralized  LRT  and  cannot  perform  as  well  as  the  two  sensor  centralized  LRT. 

In  both  these  figures  we  see  that  the  performance  initially  improves  as  A 
increases  and  then  degrades.  This  is  basically  due  to  the  effect  of  frequency 
diversity.  We  know  [ 2 J  that  most  of  the  signal  energy  is  associated  with  the 
first  2A+1  eigenfunctions  $^(t).  Thus  as  A  increases  the  number  of  signifi¬ 
cant  independent  "observations"  of  the  signal  increases,  however,  the 


*For  this  problem  the  centralized  performance  of  one  sensor  receiving  Ec  is 
exactly  equivalent  to  that  of  two  sensors  each  receiving  Ec/2. 
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signal-to-noise  ratio  (SNR)  for  each  "observation"  decreases.  Thus  these 
figures  plot  the  tradeoff  of  number  of  observations  versus  SNR  per 
observation. 

We  note  from  Figs.  3-4  and  3-5  that  the  performance  of  the  optimal  DLRT 
seems  to  be  more  similar  to  that  of  the  one  sensor  centralized  LRT  than  to 
that  of  the  two  sensor  LRT.  In  Fig.  3-6  we  consider  the  relative  performance 
of  these  systems  more  closely.  We  compare  the  centralized  and  decentralized 
systems  by  determining  the  amount  of  signal  energy  required  by  a  centralized 
system  to  perform  as  well  as  a  decentralized  system.  The  ratio  Ec/Ed  is 
largest  when  the  decentralized  system  performs  well  and  is  smallest  when  it 
performs  poorly.  We  take  as  baseline  the  performance  of  the  optimal  DLRT  with 
Ed  -100  and  Pp=0.001  and  determine  the  Ec  required  to  obtain  the  equivalent  P^. 
Note  that  for  A“2,  Ec/Ed=1.40  and  thus  the  two-sensor  DLRT  performs  as  well  as 
1.40  centralized  sensors.  As  4  increases  the  DLRT  performance  degrades  until 
at  A* 100  we  have  Ec/EdSl.25. 

The  performance  of  the  DLRT  begins  to  improve  as  2WT  increases  above  100. 
The  asymptotic  performance  has  not  been  determined  as  optimal  performance  has 
not  been  determined  analytically  and,  for  2WTM00,  numerical  difficulties  in¬ 
trude  into  the  computation  of  P^  and  Pp. 

In  Fig.  3-7  we  plot  the  performance  of  DLRTs  for  different  combinations 

i 

of  sensors  with  Pp“0.001.  If  we  let  E<j  denote  the  signal  energy  received  by 

1  2 

the  i-th  sensor  than  for  all  curves  we  have  E  +E  “200.  The  different  curves 

d  d 

1  2 

correspond  to  various  ratios  E  /E  .  We  note  that  the  performance  increases  as 

d  d 

the  sensors  become  more  asymmetric.  This  occurs  since  we  are  effectively 

1  2 

moving  toward  a  centralized  solution  (note  that  at  E  “0,  E  “200  the 

d  d 
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performance  is  identical  to  that  of  Ec»200).  When  E^O  and  E2«200  the  first 

d  d 

threshold  is  zero  so  that  u^l  always  and  the  second  threshold  is  identical 
to  that  associated  with  Ec*200. 
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SECTION  4 

SUMMARY,  CONCLUSIONS,  AND  SUGGESTED  RESEARCH 

In  this  paper  we  have  considered  the  problem  of  distributed  detection 
with  limited  communication.  We  have  obtained  optimal  detection  laws  for  the 
case  of  known  signals  in  noise  uncorrelated  between  sensors  and  for  the  case 
of  linearly  dependent  known  signals  in  white  noise  correlated  between  sensors. 
For  both  of  these  cases  the  optimal  distributed  detection  law  consists  of 
forming  local  likelihood  ratios  and  testing  these  ratios  against  thresholds  to 
determine  a  decision.  For  uncorrelated  noise  a  single  threshold  is  optimal. 
This  may  well  be  the  case  for  correlated  noise  as  well,  but  we  were  unable  to 
prove  or  disprove  this  conjecture. 

For  a  specific  example,  with  noise  correlated  between  sensors,  we  found 
that  the  distributed  detection  law  performs  worse  as  the  correlation  in¬ 
creases,  while  the  centralized  law  may  either  perform  better  or  worse.  Thus 
an  analysis  is  necessary  in  each  case  to  determine  whether  the  performance 
penalty  for  implementing  a  distributed  detection  law  is  worth  the  payoff  in 
terms  of  reduced  communications  requirements. 

We  then  investigated  the  problem  of  distributed  detection  of  an  unknown 
signal  in  noise.  This  problem  was  found  to  be  considerably  more  difficult 
than  the  case  of  a  known  signal  since  the  detection  law  depends  upon  the  en¬ 
tire  received  signal  -  a  scalar  sufficient  statistic  does  not  exist  in  gener¬ 
al.  As  the  optimal  distributed  detection  law  could  not  be  determined  we 
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investigated  laws  employing  local  likelihood  ratio  tests.  Even  for  this  case 
results  are  difficult  to  obtain  since  the  local  log-likelihood  ratios  are  not 
Gaussian.  However,  for  an  important  special  case  (long  observation  time, 
ideal  bandlimited  signals,  white  sensor  noise)  we  were  able  to  obtain  analyt¬ 
ical  results.  For  this  case  we  found  that  two  identical  distributed  sensors 
generally  perform  as  well  as  1.25  to  1.40  centralized  sensors. 

More  general  problems  of  unknown  signals  in  noise  do  not  appear  solvable 
by  using  the  methods  of  this  report  since  the  exact  probability  distribution 
function  is  generally  difficult  to  determine.  The  approach  used  in  central¬ 
ized  problems  [1),[2]  is  to  generate  approximations  based  on  Chernov  bounds 
and  Gaussian  approximations  or  Edgeworth  expansions.  We  considv _£veral 
such  approximations  for  the  distributed  problem  -  in  all  cases  the  approxima¬ 
tions  were  too  inaccurate  to  be  useful.  A  further  drawback  to  those  approxi¬ 
mations  is  that  the  moment  generating  function  m(r,s)  is  required  and,  as 
Appendix  A  makes  clear,  these  quantities  are  difficult  to  compute  for  all  but 
the  ideal  bandlimited  case. 

Work  in  distributed  detection  is  in  its  infancy,  and  many  of  the  issues 
that  arise  in  the  centralized  theory  can  also  be  formulated  in  the  distributed 
setting,  floreover,  the  increasing  tendency  to  net  sensors  together  in  mili¬ 
tary  surveillance  systems  makes  the  distributed  detection  problem  formulation 
of  potential  practical  significance.  However,  as  the  results  of  this  report 
make  clear,  the  distributed  version  of  a  given  detection  problem  is  often  sig¬ 
nificantly  more  difficult  to  solve  than  its  centralized  counter  part. 
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APPENDIX  A 

DERIVATION  OF  LRT  PROBABILITY  DISTRIBUTIONS 


In  this  appendix  we  derive  the  joint  probability  distribution  of  the  log- 
likelihood  ratios  considered  in  Eq.  3-8.  Recall  that  the  log-likelihood 
ratios  A ^  are  defined  as 

K 

Ai  •  lim  Ai  (A-l) 

K-*-* 


where 


K  1  K 
Ai  A  -  I 
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(A-2) 


and  the  are  the  eigenvalues  of  Ks(t,u)  where 

Ks( t ,u)  -  E{s(t)s(u) }  .  (A-3) 

k 

The  y±  are  the  Karhunen-Loeve  coefficients  associated  with  the  k-th  eigenfunc¬ 
tion  of  Ks.  It  is  straightforward  to  verify  that 
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where  A^  are  the  eigenvalues  associated  with  K(t,u),  then  for  long  enough  ob¬ 
servation  time  T 


gx  a  t  /  g(S(f))df 


(A-14) 


where  S(f)  is  the  spectrum  of  K(t,u).  The  can  be  written  as  in  Eq.  A-13  by 

using  Eq.  A-14  if  we  first  take  the  logarithm  and  the  limit  as  K-*».  Thus 
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where 
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We  define  Ej«c2T  as  the  signal  energy  received  by  the  i-th  sensor,  A*2WT 
as  the  observation  time-signal  bandwidth  product,  and  Yi«c2T/2WT  as  the  sig- 
nal-to-noise  ratio  in  the  signal  bandwidth.  With  these  definitions  we  have 
obtained  Eqs.  3-15  and  3-16  of  Section  3: 
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Using  a  standard  set  of  Laplace  transforms  [5]  we  can  invert  Eq.  A-21  to 
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Inverting  Eq.  A-22  on  r  first  yields 
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(continued) 
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where  Iv(»)  is  the  modified  Bessel  function  of  order  v. 

To  obtain  the  distribution  functions  we  now  integrate  Eqs.  A-23  and  A-26. 
To  perform  these  integrations  we  assume  that  A/2-1  is  a  nonnegative  integer. 

We  obtain  under  H°: 
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where 

1+Yl 

«i  -  -  Ti  ,  i-1,2  .  (A-28) 

Ti 

To  integrate  p(l1,£2|H1)  we  expand  Iv(*)  in  a  series  and  integrate  term 
by  term.  From  1 1 3 J  we  have  that 
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and  thus  we  obtain 
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(continued) 
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where 


0i 
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(A-32) 


We  thus  have  derived  closed-form  expressions  for  the  joint  probability  distri¬ 
bution  functions  as  desired. 


